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Abstract 

The Stirling refrigeration cycle using an ideal Bose-gas as the working substance is called 
the Bose-Stirling refrigeration cycle, which is different from other thermodynamic cycles such 
as the Carnot cycle, Ericsson cycle, Brayton cycle, Otto cycle, Diesel cycle and Atkinson cycle 
working with an ideal Bose gas and may be operated across the critical temperature of Bose- 
Einstein condensation of the Bose system. The performance of the cycle is investigated, based 
on the equation of state of an ideal Bose gas. The inherent regenerative losses of the cycle are 
considered and the coefficient of performance and the amount of refrigeration of the cycle are 
calculated. The results obtained here are compared with those derived from the classical Stir¬ 
ling refrigeration cycle, using an ideal gas as the working substance. The influence of quantum 
degeneracy and inherent regenerative losses on the performance of the Bose Stirling refriger¬ 
ation cycle operated in different temperature regions is discussed in detail, and consequently, 
general performance characteristics of the cycle are revealed. 
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Nomenclature 

C v heat capacity at constant volume (J) 

F(z) correction function 

g weight factor 

gi(z) Bose function 

h Planck’s constant (J s) 

k Boltzmann’s constant (J/K) 

m rest mass of a particle (kg) 

N total number of particles 

N 0 particle number in the ground state 

n number density (m -3 ) 

P pressure (Pa) 

Qi amount of heat exchanged between the working substance and the 
regenerator for a constant V H process (J) 

Q 2 amount of heat exchanged between the working substance and the 
regenerator for a constant V\_ process (J) 

Qn amount of heat exchanged between the working substance and the 
heat reservoir at temperature T u during the isothermal process 

(J) 

Q L amount of heat exchanged between the working substance and the heat 

reservoir at temperature 7) during the isothermal process (J) 

R e relative coefficient of performance 

Rq l relative refrigeration load 

r v volume ratio of two isochoric processes (J) 

S entropy of the working substance (J/K) 

T temperature of the working substance (K) 

T c condensation temperature (K) 

7’ k temperature at which the heat capacities at two constant volumes are 
equal to each other (K) 

T h temperature of the high-temperature heat-reservoir (K) 

T L temperature of the low-temperature heat-reservoir (K) 

U internal energy (J) 

V volume (m 3 ) 

V H maximum volume of the gas (m 3 ) 

V L minimum volume of the gas (m 3 ) 

v average volume that an ideal Bose particle occupies (m 3 ) 

/ i i maximum average volume of an ideal Bose particle (m 3 ) 

v L minimum average volume of an ideal Bose particle (m 3 ) 

W work input per cycle (J) 

r fugacity of the gas 

A Q regenerative loss (J) 

£ coefficient of performance 

1 mean thermal wavelength (m) 
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/< chemical potential of the gas (J) 

t temperature ratio of two heat-reservoirs 

r(I) Gamma function 


1. Introduction 

It is well known that the performance of a thermodynamic cycle using a classical 
ideal gas as the working substance can be found directly from thermodynamics text¬ 
books. However, when the gas temperature is low enough or density is high enough, 
the gas will deviate from classical gas behavior and quantum degeneracy of the gas 
will become important [1-3]. It is thus clear that one needs to consider the influence 
of quantum degeneracy on the performance of a thermodynamic cycle at low tem¬ 
peratures. Recently, several authors have analyzed the effect of quantum degeneracy 
on the performance of the engine and refrigeration cycles working with an ideal Bose 
or Fermi gas [4-9], based on quantum statistical mechanics. The influence of several 
factors on the performance of these cycles has been analyzed and many meaningful 
conclusions have been obtained. 

Like classical thermodynamic cycles [10,11], quantum thermodynamic cycles 
using an ideal Bose gas as the working substance may have different cycle models 
such as the Carnot cycle, Ericsson cycle, Stirling cycle, Brayton cycle, Otto cycle, 
Diesel cycle, Atkinson cycle, and so on. It has been proved [12], according to the 
properties of an ideal Bose gas, that some important thermodynamic processes such 
as a constant pressure or an adiabatic process cannot be carried out from the tem¬ 
perature region of T > T c to that of T < T c , where T is the temperature of the Bose 
system and T c is the critical temperature of Bose-Einstein condensation (BEC) of the 
Bose system. Consequently, the Carnot cycle, Brayton cycle, Otto cycle, Ericsson cy¬ 
cle, Diesel cycle and Atkinson cycle working with an ideal Bose gas can only be oper¬ 
ated in the temperature region above T c because these thermodynamic cycles include 
the adiabatic or/and isobaric process. The Stirling cycle is also one of the important 
thermodynamic cycles. However, unlike the thermodynamic cycles mentioned above, 
the Stirling cycle consisting of two isochoric and two isothermal processes does not 
include the adiabatic or isobaric process and can be operated across the critical tem¬ 
perature of BEC of the Bose system when an ideal Bose gas is used as the working 
substance. Thus, the Stirling cycle working with an ideal Bose gas can be operated in 
different temperature regions and will have some novel performance characteristics 
which need to be studied further. 

In the present paper, we will study the performance of the Stirling refrigeration 
cycle using an ideal Bose gas as the working substance. First, the general expressions 
for several important parameters such as the coefficient of performance, amount of 
refrigeration and work input per cycle are derived, based on the equation of state of 
an ideal Bose gas; secondly, the influence of quantum degeneracy and inherent regen¬ 
erative losses on the performance of the cycle are analyzed; Finally, some special 
cases are discussed in detail, so that the general performance characteristics of the 
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Bose Stirling refrigeration cycle operated in the various different temperature regions 
are revealed. 


2. The Bose-Stirling refrigeration-cycle model 

The Stirling refrigeration cycle working with an ideal Bose gas is called the Bose- 
Stirling refrigeration cycle, which is composed of two isothermal and two isochoric 
processes. It is operated between the heat sink at temperature T H and the cooled 
space at temperature 7\ . In order to improve the performance of the cycle, a regen¬ 
erator is used in the isochoric processes. The temperature-entropy diagram of the 
cycle is shown in Fig. 1, where Q L and Q H are the amounts of heat exchanged 
between the working substance and the heat reservoirs at temperatures T L and T H 
during the isothermal processes, Q x and Q 2 are the amounts of heat exchanged be¬ 
tween the working substance and the regenerator during the isochoric processes, and 
Ki and V H are the minimum and maximum volumes of the gas system. For the sake 
of convenience, all the heats Q L , Qh, Q\ and Q 2 are selected to be positive. 

For an ideal Bose gas, the expressions of its pressure and number density of par- 


tides are given by [1-3] 


P = gkTr 3 g 5/2 (z) 

(1) 

and 


n = (N-N 0 )/V = g r 3 g 3/2 (z), 

(2) 


respectively, where g is the number of possible spin orientations of a gas particle, k is 
the Boltzmann constant, T, N, N 0 , and V are, respectively, the gas temperature, total 



Fig. 1. The temperature-entropy diagram of a Bose-Stirling refrigeration-cycle. 
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number of particles, number of particles in the ground state, and volume, z = exp (fjJ 
kT) and / = li(2nmkT) 1/2 are, respectively, the fugacity of the gas and the mean 
thermal wavelength of particles, // is the chemical potential of the gas, h is the Planck 
constant, m is the rest mass of a particle 


gi(z) 


1 r°° x 1 - 1 dx 

W) J 0 r 


(3) 


is called the Bose function, and T(/) is the Gamma function. When T ^ T c , the par¬ 
ticle population of the ground state is macroscopically negligible compared to the 
total particle number of the ideal Bose gas system [1-3]. Consequently, it is seen from 
Eq. (2) that for an ideal Bose gas, the number density of particles, n, should be equal 
to the number density of particles that have energies larger than the energy of the 
ground state. Using Eqs. (1) and (2), one can find that the internal energy and entro¬ 
py of an ideal Bose gas system are, respectively, given by [1] 



|NkTg 5/ 2 (z)/g 3/2 (z), T Ss T c , 

\{N~N 0 )kT §fg = |p§Mr(£) 3/2 , T < T c , 


(4) 


and 


S = 


' Nk[lg 5/2 (z)/g 3/2 (z) - In z\, T ^ T c 

i(^-Wf|| = !§gMr(x) 3/2 , t <t c , 

where ((5/2) = 1.341, {(3/2) = 2.612, T c = 

From E 
pressed as 


(5) 


and w is the s P ecific volume - 

From Eq. (4), one can prove that the heat capacity at constant volume can be ex- 


Cv=(^) =^Nk^[TF(T,v)} = 




AT 1 


£5/26) _ 3 gyi(~) 

gl/l{z) 5 g 1/2 to 

i5 M «V3)fxV V2 

4 lyK £(3/2) \t c J > 


T ^ T c 


T < T c 


( 6 ) 


where 


F(T,v) = 


( gs/ll-CT, v)\/g 3/2 M T , ”)]> T > T e , 

1 [C(5/2)/C(3/2)][7’/7’ c ] 3/2 , T<T c . 


Using Eqs. (5) and (6), the amounts of heat exchanged in an isothermal and an 
isochoric process are, respectively, given by 


ef f = 


[ 1 TdS = T(S{ - S{) 

JSi 


= NkT{ - [F(T, Vi) - F(T, Hi)] - [In z(T, v { ) - In z(T, r,)[ 


(V) 
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and 

a f = J T ' c v (r, ®) dr = ^M[r f F(r f , w ) - »)], (8) 

where the indices i and f refer to the initial and final states, respectively. 


3. Regenerative characteristics 

From Eq. (8), one can find that the net amount of heat transfer between the work¬ 
ing substance and the regenerator during two isochoric processes is determined by 

A Q=Q 1 ~Q l = [ H C v (7> L )dr- f H Cy(T, r H )dr 

= ^M{r H [E(r H ,t; L ) -F{T H ,V H )} -T l [F{T l ,v l )-F(T l ,v h )}}. (9) 

Using Eq. (6), we can plot the CylNk versus T characteristic curves, as shown in 
Fig. 2. It is clearly seen from Fig. 2 that Cy(T, u L ) > Cy(T,v H ) when T> T E and 
Cy(T, dh) > Cy(T,v L ) when T < T E , where T E is the temperature at which the heat 
capacities at two constant volumes are equal to each other, i.e., Cy(T E , v E )- 
C v (T e , l' h ). Consequently, it can be seen from Eq. (9) that there are two possible 
cases: (a) A Q > 0 and (b) A Q < 0. When A Q > 0, the amount of heat (9 2 flowing into 
the regenerator in the small constant-volume regenerative process is larger than that 
of Q i flowing from the regenerator in the large constant-volume regenerative pro¬ 
cess. The redundant heat in the regenerator per cycle must be released to the cold 
reservoir in a timely manner. This results in the reduction of the amount of refriger¬ 
ation from Q e to 0 l- ^ not, the temperature of the regenerator would be changed 



T(K) 


Fig. 2. The curves of Cy varying with T for two different specific volumes. 
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such that the regenerator would not be operated normally. Similarly, when AQ < 0, 
the amount of heat Q 2 flowing into the regenerator in the small constant-volume 
regenerative process is smaller than that of Q\ flowing from the regenerator in the 
large constant-volume regenerative process. The inadequate heat in the regenerator 
per cycle must be compensated from the hot reservoir in a timely manner, while the 
amount of refrigeration Q L is unvarying. 

According to the regenerative characteristics mentioned above, the unified expres¬ 
sion for the amount of refrigeration per cycle can be given by 

& = Ql~ <5A Q 

= NkT^-[F{T u Vn) -F{T u vQ)\ - [lnz(T L , u H ) - lnz(r L , «l)]| 

- S^Nk{T H [F(T Hl v L ) -F(r H ,PH)] - T l [F(T l ,v l )-F(T l ,v h )}}, (10) 
where d - 1 when AQ > 0 and <5 = 0 when AQ < 0. 

4. Expressions for several important parameters 

Besides the amount of refrigeration, the coefficient of performance and work in¬ 
put are also the important performance parameters of the cycle. Using Eqs. (7), (8) 
and (10), one can find that the work input per cycle and coefficient of performance 
may be, respectively, expressed as 

W = Q h - 0L + 02 - Ql 

= NkT\\[F(T H , Uh) — F(Th, ^l) + lnz(7’ H , i-’l) — lnz(7’ H , Uh)] 

+ NkT H [F(T H , v H ) - F(T h ,v l ) + In z(T H , v L ) - Inz(r H , v H )] (11) 

and 

T L {f [E(z d ) - F(z c )] - ln(z d /z c )} - 8\{T n [F{z h ) - F(z a )] - T L [F(z c ) - F(z d )]} 
Th[F(zQ) - F(z b ) + ln(z b /z a )] + T L [F(z c ) - F(z d ) + ln(z d /z c )] 

(12) 

where z a , z b , z c and z d are the fugacity of the gas in state points a, b, c and d shown in 
Fig. 1. 

In addition, in order to compare the performance of a Bose-Stirling refrigeration 
cycle with that of a classical Stirling refrigeration cycle using a classical ideal gas as 
the working substance, we introduce the relative amount of refrigeration and relative 
coefficient of performance as follows: 


| [F(z d ) - F(z c )\ - ln(z d /z c ) - d|{T[E(z b ) - F(z a )] - [F{z c ) - F(z d )]| 

lnr v 


( 13 ) 
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and 


-£ 

s c 

_ (r - l){|[F(z d ) -F(z c )\ - In(z d /z c ) - 5\{x [F(z b ) - F(z a )] +F(z d ) - F(z c )}} 
t[F(z a ) - F(z b ) + ln(z b /z a )] + F(z c ) - F(z d ) + ln(z d /z c ) 


(14) 

where Q\ = AT: 7) In r v and e c = 7j/( F H — 7' [ ) are, respectively, the amount of refrig¬ 
eration and coefficient of performance of a classical Stirling refrigeration cycle with 
perfect regeneration, and x = 7 h /7'i. and r v = are, respectively, the temperature 
ratio of the reservoirs and the volume ratio of two constant-volume processes. 

Using the above equations, we can discuss the performance characteristics of the 
cycle operated in the various different temperature-regions. 


5. General performance characteristics 

It is clearly seen from Eqs. (9) (12) that, for a Bose-Stirling refrigeration cycle, the 
amount of refrigeration, work input and coefficient of performance are, in general, 
dependent on the regeneration characteristic, temperature, volume and other param¬ 
eters, while the net heat addition to or rejection from the regenerator depends on the 
magnitude order of T H , 7j and T E except for depending on the temperature and vol¬ 
ume. It is obvious that there are different performance characteristics for the Bose- 
Stirling refrigeration cycle operated in the different temperature regions, which will 
be discussed below. 


• When T h > T l > T b , Cy(T, r H ) < C\-(T. v L ) and AQ > 0. For the case of 
T l > 7’ c (/.'| ), Eqs. (12) (14) can be, respectively, written as 

... 7\[F(z d ) - F(z c ) - ln(z d /z c )] + (3/2)?ffi[F(z a ) - F(z b )\ 

T\i[F(z a ) - F(z b ) + ln(z b /z a )] - 7 n L [F(z d ) - F(z c ) - ln(z d /z c )] ’ 

R _ F{z a ) - F(z c ) - ln(z d /z c ) + (3/2)t[F(z a ) - F(z b )] 

QL In r v 

and 

R _ (t - l){F(z d ) - F(z c ) - ln(z d /z c ) + (3/2)r[F(z a ) - F(z b )]| 
r[F(z a ) - F(z b ) + ln(z b /z a )] + F(z c ) - F(z d ) + ln(z d /z c ) 

while for the case of T c (vy) > 7j, Eqs. (12)—(14) can be further simplified as 

7\[F(z d ) - 0.5134(7Vr CL ) 3/2 - lnz d ] + (3/2)r H [F(z a ) - F(z b )] 

Tn[F{z & ) - F(z b ) + ln(z b /z a )] - T L [F(z d ) - 0.5134 (T L /T ch ) y2 - lnz d ] ’ 


(15) 

(16) 

(17) 


F(z d ) - 0.5134(7Y/r CL ) 3/2 ~ lnz d + (3/2 )t[F(z a ) - F(z b )] 


In r, 


(18) 

(19) 
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and 


(t - l){F(z d ) - 0.5134(r L /r CL ) 3/2 - Inz d + (3/2)r[F(z a ) - F(z b )]} 
x[F(z a ) - F(z h ) + ln(z b /z a )] + 0.5134(7’ L /7’ CL ) 3/2 --F(z d ) + ln(z d ) 

( 20 ) 


where Tcl = 


27lmA'[i' L 5(3/2)] 2/3 ‘ 


Using Eqs. (1)—(3), (16) and (19) and choosing He 4 gas as the ideal Bose gas, one 
can plot the curves of R QL varying with the temperature 7) of the cooled space for 
given u L , t and r v , as shown in Fig. 3. It is clearly seen from Fig. 3 that when the 
temperature ratio of the heat reservoirs t is small, the relative amount of refrigera¬ 
tion Rq\ is a monotonically decreasing function of T L , but it is always greater than 
unity; when the temperature ratio of the heat reservoirs x is large, 7 ?q L is a monoton¬ 
ically increasing function of 7) , but it is always smaller than unity. This implies that 
when t is small, the amount of refrigeration of the Bose-Stirling refrigeration cycle is 
greater than that of the classical Stirling refrigeration cycle because the effect of 
quantum degeneracy on the amount of refrigeration is larger than that of regenera¬ 
tive losses; when x is large, the amount of refrigeration of the Bose-Stirling refriger¬ 
ation cycle is smaller than that of the classical Stirling refrigeration cycle because the 
effect of regenerative losses on the amount of refrigeration is larger than that of 
quantum degeneracy. In addition, Fig. 3 shows that there is a value of the temper¬ 
ature ratio t for which the amount of refrigeration ratio 7? QL is equal to unity. This 
indicates that the influence of quantum degeneracy and regenerative losses on the 
amount of refrigeration of the cycle counteract one another in this case, so that 
the amount of refrigeration of the Bose-Stirling refrigeration cycle is the same as that 
of the classical Stirling refrigeration cycle. It can be also seen, from Fig. 3, that when 
the temperature ratio t is small, the larger the volume ratio r v is, the smaller the rel- 



T l (K) 


Fig. 3. The curves of Rql varying with T L in the temperature region of T H > T L > T E for some given 
values of v E , r v and z. 











522 


B. Lin et al. / Applied Energy 83 (2006) 513-535 


ative amount of refrigeration i? QL ; when the temperature ratio x is large, the larger 
the volume ratio r v and the larger the relative amount of refrigeration i? QL . Obvi¬ 
ously, R Qh tends to unity with the increase of T L . This is an expected result because 
the quantum behavior of gas particles is negligible in this case. 

Similarly, using Eqs. (1) (3), (17) and (20) and choosing He 4 gas as the ideal Bose 
gas, we can plot the curves of R e varying with the temperature 7) of the cooled space 
for given different parameter conditions, as shown in Fig. 4. It can be seen, from 
Fig. 4, that the relative coefficient of performance R, is a monotonically increasing 
function of T { for given v L , r and r v , but it is always smaller than unity. This implies 
that the coefficient of performance of the Bose-Stirling refrigeration cycle is smaller 
than that of the classical Stirling refrigeration cycle and decreases as the quantum 
degeneracy of the working substance increases. The figure shows that for given 
T l , v i and r v , the larger the temperature ratio x, the smaller the relative coefficient 
of performance R e ; for given T L , t> L and x, the larger the volume ratio r v , the larger 
the R,. 

• When T h > T c (v t ) > T v > 7) > 7’ c (r n ). there may have two cases of A Q > 0 and 
A Q < 0. If/; r "(C V L - C VH )dr > f£(C VH - C VL )dr, AQ > 0 and Eqs. (12) (14) can 
be, respectively, expressed as 


T l [F{z a ) - 0.5134(T L /r CL ) 3/2 - lnz d ] + (3/2 )T n [F(z a ) -F(= h )} 
T H [F(z a ) - F(z b ) + ln(z b /z a )] - T L [F(z d ) - 0.5134(r L /r CL ) 3/2 - lnz d ] 

_ F{z & ) - 0.5134(T l /7 , C l) 3/2 - In z d + (3/2)r[F(z a ) - F(z b )] 

In (uh/ul) 


( 21 ) 

( 22 ) 


and 



Fig. 4. The curves of R e varying with J L in the temperature region of 7 h > T l > T e for some given values 
of v L , r v and t. 
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(t - l){F(z d ) - 0.5134(r L /r CL ) 3/2 - Inz d + (3/2)r[F(z a ) - F(z b )]} 

r[F(z a ) - F(z b ) + ln(z b /z a )] + 0.5134(r L /r C L) 3/2 - F(z & ) + ln(z d ) 


(23) 

If / r 7 E H (C VL -CvH)dr < f(((CyH — Cvl)&T, /±Q<0 and Eqs. (12) (14) can be, 
respectively, expressed as 

r L {(5/2)[F(z d ) - 0.5134(r L /r CL ) 3/2 ] - lnz d } 

T H [F(z a ) - F(z b ) + ln(z b /z a )] - T L [F(z A ) - 0.5l34(T L /T CL f 2 - lnz d ] ’ 

(24) 


(5/2 )[E(z d ) - 0.5134(r L /r CL ) 3/2 ] - lnz d 
ln(r H /r L ) 


(25) 




Fig. 5. The curves of Rq E varying with r v in the temperature region of T H > T c (i’l) > T e > T e > T c (vji) for 
some given values of v E , T e and r. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 
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and 

(t - l){(5/2)[F(z d ) - 0.5134(:T l /Tcl) 3/2 ] - lnz d } 
t[F(z a ) - F(z b ) + In (z b /z a )] - [F(z d ) - 0.5134(F l / F C l) 3/2 - lnz d ]' 


(26) 


When AQ > 0, the variations of the relative amount of refrigeration Rq\ and the 
relative coefficient of performance R £ with the volume ratio r v for given values of 7). 
v L and t are shown in Figs. 5(a) and 6(a), respectively. It is seen, from the figures, 
that for given 7\ and r L , when the temperature ratio r is small, R<,\ is larger than 
unity and it is a monotonically decreasing function of r v , while the coefficient of per¬ 
formance of the Bose-Stirling refrigeration cycle is smaller than that of the classical 
Stirling refrigeration cycle and it is a monotonically increasing function of /%.; when 
the temperature ratio t is large, both R ()l and R, are always smaller than unity and 
both of them are monotonically increasing functions of i\. For given T L , v L and r v , 
the larger the temperature ratio i, the smaller the Rq E and R y . 

When AQ < 0, we can obtain the curves of both Rq L and R. varying with the vol¬ 
ume ratio /\, for given different parameter conditions, as shown in Figs. 5(b) and 6(b), 
respectively. It is seen from the curves in the figures that both Rq E and R e are smaller 
than unity and they are monotonically increasing functions of r v for given 7) , v { and 
t. This indicates that the amount of refrigeration and coefficient of performance of 
the Bose Stirling refrigeration cycle are smaller than those of the classical Stirling 
refrigeration cycle. 

• When T u > T c (v L ) > T E > T c (v H ) > 7i , there may have two cases of AQ > 0 and 
AQ < 0. If fr E H (Cy L - C VH ) d T > j£(C v H - C VL ) d T, AQ>0 and Eqs. (12)-(14) are, 
respectively, expressed as 


o.5i34r L [(r L /r CH ) 3/2 - (t l /t cl ) 3/2 } + ( 3 / 2 )r H [F(z a ) -F(r b )] 
Tn[F(z a ) - F(z b ) + ln(zb/z a )] - 0.5134r L [(T L /r CH ) 3/2 - (F L /r CL ) 3/2 ] 

_ Q.5i34[(r L /r CH ) 3/2 - (r L /r CL ) 3/2 ] + (3/2)t[F(z a ) - F(z b )] 

ln(u H /fL) 


(27) 

(28) 


and 

(t - l){0.5134[(F L /r CH ) 3/2 ~ (Fl/Fcl) 372 ] + (3/2)r[F(z a ) - F(z b )]} 
t[F(z a ) - F(z b ) + ln(z b /z a )] - 0.5134[(7 ’ l /F ch ) 3/2 - (T h /T CE f 2 } 

(29) 

Similarly, if /^“(Cvl — C V h) dF < /^ r (C V H — C V l) dF, AQ < 0 and Eqs. ( 12) (14) 
can be, respectively, expressed as 


_1.2835F L [(r L /F CH ) 3/2 - (T L /T CL f 2 ] _ 

F H [F(z a ) - F(z b ) + ln(z b /z a )] - 0.5134r L [(T L /r CH ) 3/2 - (F L /r CL ) 3/2 ] 


1.2835[(r L /Tc H ) 3/2 -(r L /FcL) 3/2 ] 

ln(v H /v L ) 


(30) 

(31) 
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Fig. 6. The curves of R e varying with r v in the temperature region of Ta > T c (!; l ) > T e > T l > T c (v h ) for 
some given values of v E , T h and r. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 


and 

R= _ 1.2835(t - i)[(r/r CH ) 3/2 - (TJT CL ) 3/2 ] _ 

t[F(r a ) - F(z b ) + ln(r b /z a )] - 0.5\34[(T L /T CH f 2 - (T L /T CL f 2 } ’ 
where T cH = - 

Using above equations, one can plot the curves of f? QL and R t . varying with 
the volume ratio r v for given r L , T { and t, as shown in Figs. 7 and 8, respec¬ 
tively. It is seen from the curves in Figs. 7(a) and 8(a) that, in order to satisfy 
the condition of A Q > 0, the temperature ratio t has to be chosen to be very large 
and that for given u L , 7) and x, both /? QI and II, are monotonically decreasing 
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Fig. 7. The curves of Rq L varying with r v in the temperature region of Th > T c (v L ) > T e > T c (v h ) > T l for 
some given values of r L , T L and t. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 


functions of r v and the larger the temperature ratio t, the smaller the Rq\ (and 
R,). This implies that the effect of the regenerative losses on the performance of 
the cycle is much larger than that of quantum degeneracy and, in such a case, the 
larger the volume ratio r v , the larger the regenerative losses. Figs. 7(b) and 8(b) 
show that when AQ < 0, both Rq L and R e are monotonically increasing functions 
of r v for given v L , 7) and r. This indicates that when A Q < 0, the effect of quan¬ 
tum degeneracy on the performance of the cycle is much larger than that of 
regenerative losses and the larger the volume ratio r v , the stronger the quantum 
degeneracy. 

• When T c (v L ) > T H > T E > T L > T c (v H ), there are two cases AQ > 0 and AQ < 0. 
When fr“(C VL - C VH )dr > f^(C VH - C VL )dT, AQ > 0 and Eqs. (12)-(14) are, 
respectively, written as 
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Fig. 8. The curves of R t: varying with r v in the temperature region of T H > TJpiO > T E > T c (v h) > Ty for 
some given values of v L , T L and r. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 


S = 


Ty 

/ \ 3/2 

F(z d )-0.5134^J -lnz d 

+ |r H 

F{ z a) — 0.5134 ( 



Th 

/ \ 3 / 2 

F(z a )-0.5134(^-J — lnz a 

- T L 

r(z a )-0.5134(i.; 

3/2 

) - lnz d 


(33) 


F(z d ) - 0.5134 (T l /T cl ) 3/1 - lnz d + (3/2)t[F(z a ) - 0.5134(r H /r CL ) 3/2 ] 

ln(r H /t)L) 


( 34 ) 





528 


B. Lin et al. / Applied Energy 83 (2006) 513-535 


and 


(t - l){F(z d ) - 0.5134(7 L /7 CL ) 3/2 - lnz d + 1.5t[F(z a ) - 0.5134(r H /r CL ) 3/2 ]} 

i[F(z a ) - 0.5134(7 H /7 CL ) 3/2 - lnz a ] +0.5134(r L /r CL ) 3/2 -F(z d ) + ln(z d ) 

(35) 

While /r E H (C V L - C VH )d7 < J^(C VH - C VL )d7, AQ < 0 and Eqs. (12)-(14) can be 
further simplified as 

7 L {(5/2)[F(z d ) - 0.5134(7’ l / 7’ cl ) 3/2 ] - lnz d } 

T*[F{z a ) - 0.5134(7’ h / 7’ cl ) 3/2 - lnzj - 7 L [F(z d ) - 0.5134(7 l /7 C l) 3/2 - lnz d ] ’ 

(36) 


Rqy = 


R,= 


(5/2)[F(z d ) - 0.5134(r L /r C L) n - lnzd 
ln(r H /rL) 

(t- l){(5/2)[F(z d )-0.5134(r L /r CL ) 3/2 ]-lnz d } 


r[F(z a ) - 0.5134(r H /r CL ) 3/2 - lnzj - [F(z d ) - 0.5134(7 L /7 CL ) 3/2 - lnz d ] 


3/2 


(37) 


(38) 


In this case, the variations of Rq E and R s , with r v are shown in Figs. 9 and 10. It is 
clearly seen from Figs. 9 and 10 that, for the case of AQ > 0, R e is always smaller 
than unity, while Rq E is always larger than unity and both R QL and R. are monoton- 
ically decreasing functions of r v ; while for AQ < 0, both R QE and R. : are monotoni- 
cally increasing functions of r v . It shows that, in this temperature order, the effect of 
quantum degeneracy of the working substance is distinct and the regenerative losses 
increase with the volume ratio r v . 

• When T c (v L ) > T H > T E > TJv H ) > 7), it is clearly seen from Fig. 2 that 
fj " (C V L - C VH ) dF < // l e (CVh — Cvl) dr. In this case, AQ < 0 and Eqs. (12) (14) 
can be, respectively, expressed as 


_i.2835r L [(r L /r CH ) 3/2 - (7 L /r CL ) 3/2 ]_ 

F H [F(z a ) - 0.5i34(r H /r CL ) 3/2 - lnzj - o.51347 l [(7 l /7 C h) 3/2 - (Tl/Tcl) V2 } ’ 
1.2835 [(7 L /r CH ) 3/2 - (F L /r CL ) 3/2 l 
ln(t> H /r L ) 

(39) 


and 


R = _ 1.2835(t - l)[(F L /r CH ) 3/2 - (F L /r CL ) 3/2 ] _ 

t[F(z a ) - 0.5134(7 ’ h /7’cl) 3/2 - lnzj - 0.5134[(7 L /7 CH ) 3/2 - (T L /T CL f 2 ]' 

(40) 

• When 7 k > 7/ r > T E > TJv H ) or T v > 7 H > 7 c (r H ) > 7), it can be seen, from 
Fig. 2, that C Y (T,v H ) < C v (T,v L ), so AQ < 0. When T E > T u > T E > T c (v n ), Eqs. 

(12)—(14) can be, respectively, expressed as 
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Fig. 9. The curves of Rq l varying with r v in the temperature region of T c (v L ) > T n > T E > T L > T c (vjj) for 
some given values of u L , T L and r. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 


r L {(5/2)[F(z d ) - 0.5134(r L /r CL ) 3/2 ] - lnz d | 

Th[F(z 3 ) - 0.5134(r H /r CL ) 3/2 - lnzj - r L [F(z d ) - 0.5134(r L /r CL ) 3/2 


(5/2)[F(z d ) - 0.5134(r L /r C L) V2 ] - lnz d 
In (v H /v L ) 


lnz d ] ’ 

(41) 

(42) 


and 


(t- l){(5/2)[F(z d ) -0.5134(r L /r CL ) 3/2 ] -lnz d } 

t[F(z a ) - 0.5134(r H /r CL ) 3/2 - lnz a ] - [F(z d ) - 0.5134(r L /r CL ) 3/2 - lnz d ]' 

(43) 
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Fig. 10. The curves of R r varying with r v in the temperature region of T c (v L ) > Th > T e > J L > r c (r H ) for 
some given values of ul, T l and x. Curves are presented for (a) A Q > 0 and (b) A Q < 0. 


When T E > T H > T c (v H ) > 7), Eqs. (12) (14) are, respectively, expressed as 

1.2835r L [(T L /r CH ) 3/2 - {T l /T C l) V2 } 


Rql — 


Th[f(z,) - 0.5134(r H /r CL ) 3/2 - inz a ] - o.5i34r L [(:r L /W /2 - (t l /t cl ) 3/2 } 
1.2835 [(r L /r CH ) 3/2 - {t l /t cl ) 3/2 } 


and 


R : = 


ln(u H /uL) 

1 .2835(t -1) [(r L /r CH ) 3/2 - (7Y/W /2 ] 


(44) 

(45) 


t[F(z a ) - 0.5134(r H /r CL ) 3/2 - lnz a ] - o.5i34[(r L /r CH ) 3/2 - (T L /T Ch f 2 }' 

(46) 

















B. Lin et al. / Applied Energy 83 (2006) 513-535 


531 


In the last two cases, the variations of i? QL and R £ with r v are shown in Figs. 11- 
16, respectively. It can be seen, from Figs. 11-16, that both R QL and R, are mono- 
tonically increasing functions of r v and their values are smaller than unity. This indi¬ 
cates that in these temperature regions, the coefficient of performance and amount of 
refrigeration of the cycle are smaller than those of the classical Stirling refrigeration 
cycle and the regenerative losses and amount of refrigeration increase with the vol¬ 
ume ratio r v . 



Fig. 11. The curves of Rq l varying with r v in the temperature region of T c (v L ) > t h >t e > t c (v h) > r L 
for some given values of v E , T E and r. 



Fig. 12. The curves of R r varying with r v in the temperature region of r c (u L ) > T H > T E > r c (f>H) > T\ for 
some given values of u L , T h and r. 
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Fig. 13. The curves of ,Rq l varying with i\ in the temperature region of t e >t h >t e > t c (k h ) f° r some 
given values of v L , T E and r. 



Fig. 14. The curves of R e varying with r v in the temperature region of t e >t h >t l > t c (r H ) for some 
given values of v E , T E and t. 


• When T c (v c ) > T E > T c (v H ) > T H > T L , it is clearly seen from Fig. 2 that 
Cy(T, i H ) > Cy(T,v L ), so A Q < 0. In this case, Eqs. (12)—(14) can be simplified as 


£ = 


5/2 

t 5/2 _ l ’ 


(5/2)[(F L /r C H) 3/2 -(F L /rcL) 3/2 ] 

In (vh/vl) 


5(r-l) 

2(r 5 /2- !)■ 


(47) 

(48) 


(49) 
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Fig. 15. The curves of Rq l varying with r v in the temperature region of T e >T h >T c (k h ) > r L for some 
given values of v E , If and t. 



Fig. 16. The curves of R e varying with r v in the temperature region of T E > T H > T c (vh) > J L for some 
given values of u L , T h and t. 


• When the temperature of the gas is high enough and its density is low enough, 
the fugacity of the Bose gas z is much smaller than unity. In such a case, g{z) = z, 
F(z) = 1, and an ideal Bose gas becomes an ideal classical gas. Eqs. (9), (10), (12), 
(13) and (14) can be, respectively, simplified as A Q - 0, Q' L = Q\ , e = s c , R y = 1 
and Rql = 1. The results are just those of a classical Stirling refrigeration cycle. This 
shows clearly that, at high temperatures, the Bose-Stirling refrigeration cycle be¬ 
comes the classical Stirling refrigeration cycle. 
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6. Conclusions 

It is clearly seen from the results obtained above that because there are inher¬ 
ent regenerative losses, the Stirling refrigeration cycle working with an ideal Bose 
gas may not possess the condition of perfect regeneration. On the other hand, 
due to the existence of BEC, the performance of the Stirling refrigeration cycle 
operated in the different temperature regions is obviously unique. The influences 
of the inherent regenerative losses and the quantum degeneracy on the coefficient 
of performance and amount of refrigeration of the cycle have been analyzed in 
detail. The coefficients of performance of the Bose-Stirling refrigeration cycle, un¬ 
der different conditions, are derived analytically. In general, the coefficient of per¬ 
formance depends not only on the temperatures of two heat-reservoirs but also 
on the volumes of two constant-volume processes and other parameters, and 
are always lower than that of a Stirling refrigeration cycle working with a classi¬ 
cal ideal gas. Only if the gas temperature is low or high enough, can the coeffi¬ 
cient of performance be a function of temperature, but independent of other 
parameters. The results obtained here reveal the general performance characteris¬ 
tics of the Bose-Stirling refrigeration cycle operated in the various different tem¬ 
perature regions. 
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